Electronic and Magnetic Structures of Chain Structured Iron Selenide Compounds 



T— I 

o 

EL, 



u 

■)—> 

a 

fl 
O 

o 



> 

O 

<N 
O 
<N 



Wei Li, 1 Chandan Setty, 2 X. H. Chen, 3 and Jiangping Hu 4 ' 2 Q 

1 Department of Physics, Fudan University, Shanghai 200433, China 
2 Department of Physics, Purdue University, West Lafayette, Indiana 41901, USA 
3 Hefei National Laboratory for Physical Science at Microscale and Department of Physics, 
University of Science and Technology of China, Hefei, Anhui 230026, China 
^Beijing National Laboratory for Condensed Matter Physics, 
and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
(Dated: February 20, 2012) 

Electronic and magnetic structures of iron selenide compounds Ce202FeSe2 (2212*) and 
BaFe2Se3(123*) are studied by the first-principles calculations. We find that while all these com- 
pounds are composed of one-dimensional (ID) Fe chain (or ladder) structures, their electronic struc- 
tures are not close to be quasi- ID. The magnetic exchange couplings between two nearest-neighbor 
(NN) chains in 2212* and between two NN two-leg-ladders in 123* are both antiferromagnetic 
(AFM), which is consistent with the presence of significant third NN AFM coupling, a common 
feature shared in other iron-chalcogenides, FeTe (11*) and K y Fe2-xSe2 (122*). In magnetic ground 
states, each Fe chain of 2212* is ferromagnetic and each two-leg ladder of 123* form a block-AFM 
structure. We suggest that all magnetic structures in iron-selenide compounds can be unified into 
an extended J1-J2-J3 model. Spin- wave excitations of the model are calculated and can be tested 
by future experiments on these two systems. 

PACS numbers: 74.25.Jb, 74.70.-b, 74.25.Ha, 71.20.-b 



I. INTRODUCTION 

The newly discovered 122*, J A 2/ Fe2- £C Se2, iron- 
chalcogenide superconductors^"— have attracted enor- 
mously interests. Like iron-pnictide high temperature 
superconductors, the FeSe-based superconductors have 
the same robust tetrahedral layers structure. How- 
ever, there are a number of distinct intriguing physical 
properties which are noticeably absent in iron-pnictide 
materials. Such as antiferromagnetically (AFM) or- 
dered insulating phases^, extremely high Neel transition 
temperatures^, and the presence of intrinsic Fe vacan- 
cies and ordering 9 12 . In addition, the Fermi surface (FS) 
topologies of superconducting compounds are very differ- 
ent from previously known superconducting Fe-pnictides. 
Both band structures calculations^"— and angle resolved 
photoemission spectroscopy studie s 17 ! 18 indicated that 
only the electron pockets are present in the supercon- 
ducting compounds, while the hole pockets around T 
point observed in iron-pncitide counterparts sink below 
the Fermi level, indicating that the inter-pocket scatter- 
ing between the hole and electron pockets is not an es- 
sential ingredient for superconductivity. 

Following the discovery of the 122* iron-chalcogenide, 
two new materials, Ce20 2 FeSe 2 (2212*)i£ and 
BaFe 2 Se 3 (123*)22r— , have been synthesized. In 
the 2212*, each Fe-layer is composed of coupled one- 
dimensional (ID) Fe-chains and in the 123*, it is 
structured by coupled two- leg ladders. An analogy of 
these materials can be made to those in cuprates, such 
as a spin-ladder system Sr^-^Ca^C^^^— »2& (x=11.5- 
15.5) and a double chain system Pr2Ba4Cu70i5_#2£. 
123*—"—, which has been investigated intensively. It was 
also reported that the 123* may be superconducting 20 . 
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FIG. 1: (Color online) Calculated crystal structures of the 
iron selenide compounds: (a) 2212*, which consists of one 
chain of FeSe4 tetrahedra structure, and (b) 123*, which con- 
sists of two- leg- ladders of edge-sharing FeSe4 tetrahedra struc- 
ture, respectively. 



In this paper, we present the theoretical study of the 
electronic band structures and magnetic orders in these 
iron selenide systems featured with low-dimensional iron 
structures. We investigate two materials including 2212* 
and 123* and show that while all these compounds are 
composed of ID Fe chain structures, their electronic 
structures are not close to be quasi- ID. Their FS still ex- 
hibit two dimensional or even three dimensional topolo- 
gies. We calculate their magnetic ordered ground states. 
In 2212*, the magnetic order ground state is a collinear- 
AFM (CAF), similar to iron-pncitides. In 123*, the 
magnetic structure is a block-AFM (BAF), similar to 
KFe2Se2 2T . The magnetic exchange couplings between 
two nearest neighbor (NN) chains in 2212* and between 
two NN two-leg-ladders in 123* are both AFM, which is 
consistent with the presence of significant third NN AFM 
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TABLE I: Structural parameters, density of states at the 
Fermi level N(Ef) (in the (eV) _1 units per Fe atom) and 
the calculated specific heat coefficient jo[mJ/(K 2 mol)], and 
Pauli susceptibility xo ( 10 ~ 9 m 3 /mol) for iron selenide com- 
pounds in NM state. The lattice parameters and the internal 
coordinates are all optimized within energy minimization. 





a(A) 


b(A) 


c(A) 


N(E F ) 


7o 


Xo 


2212* 


5.5508 


5.6794 


16.2566 


6.7293 


31.7196 


5.4546 


123* 


5.3821 


9.1123 


11.2096 


0.7650 


3.6060 


0.6201 



coupling, J3, in FeTe (11*) and 122*. This result suggests 
that all magnetic structures in iron selenide compounds 
can be unified into an extended J1-J2-J3 model. We also 
calculate spin-wave excitations of the model which can 
be tested in future experiments on these two systems. 



II. THEORIES AND RESULTS 
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FIG. 2: (Color online) Electronic structures of the iron se- 
lenide compounds in the NM state: The band structure of 
2212* (a), 123* (b), and the corresponding FS (c) and (d), 
respectively. The Fermi energy is set to zero. 



A. First-Principles Calculations 

We perform the first-principles calculations on the iron 
selenide compounds: 2212*, which has ID chains of edge- 
shared FeSe4 tetrahedra structure and 123*, which con- 
sists double chains (two-legged ladders) of edge-shared 
FeSe4 tetrahedra structure. The crystal structures are 
shown in Fig. [TJ In our calculations the projected aug- 
mented wave method 28 as implemented in the VASP 
code^ 9 -, and the Perdew-Burke-Ernzerhof exchange cor- 
relation potential 30 was used. All atomic positions and 
the lattice constants are allowed to relax simultaneously 
to minimize the energy only for nonmagnetic (NM) state. 
The experimental crystal structure a 19 i 21 are used for cal- 
culating magnetic states. A 500eV cutoff in the plane 
wave expansion ensures the calculations converge to 10 ~ 5 
eV, and all atomic positions and the lattice constants 
were optimized until the largest force on any atom was 
0.005eV/A. To properly describe the strong electron cor- 
relation in the 4/ rare earth element Ce, the LDA plus 
on-site repulsion U method (LDA+/7) was employed with 
the effective U value (U eff = U - J) of 12.0eV for the 
compound 2212*, where the U e ff value has been reported 
in the previous work 3 -^ of CeOFeAs. The results are also 
checked for consistency with varying U e ff values. We do 
not apply U e ff to the itinerant Fe-3d states. 

First, we focus on the electronic structures of the iron 
selenide compounds and their dependence on the struc- 
tural factors. For this purpose, full structural optimiza- 
tion of the these compounds were performed both over 
the lattice parameters and the atomic positions includ- 
ing the internal coordinate z of Se atom by the energy 
minimization. All these results (not included the inter- 
nal coordinate) are summarized in the Table HI Actually, 
both iron selenide compounds of the lattice parameters 
optimized in our NM calculations are found smaller by 
about 2% than the ones in experimental values^"—. In 



addition, the density of states N{Ep) at the Fermi en- 
ergy are also calculated, and the corresponding electronic 
specific heat coefficient 70 and Pauli susceptibility xo are 
all listed in Table U 

Figure [2] shows the NM state band structure and FS 
of both iron selenide compounds. As we can see that 
there are three bands crossing the Fermi level for both 
iron selenide compounds. Although both iron selenide 
compounds are composed of ID Fe chain (ladder) struc- 
tures and exhibit quasi- ID characters, their FSs still ex- 
hibit two-dimensional or even three dimensional complex 
topologies. Their NM state electronic structures are very 
distinct from that of the iron selenide superconductor 
KFe2Se2^~— . Therefore, if superconductivity exists in 
these compounds, it provides a new playground to test 
theoretical mechanisms. 

Because the NM state is strongly unstable against mo- 
ment formation, we turn to study the magnetic struc- 
tures in both iron selenide compounds. The six differ- 
ent possible magnetic configurations, as shown in Fig. [3] 
[ferromagnetic (FM) state has not been included] are all 
calculated. In Table HU we list the energies of different 
magnetic states. For 2212*, it is shown that the CAF is 
the lowest energy state. In the CAF state of 2212*, spins 
are FM for each Fe-chains and are AFM between two NN 
chains. The calculated magnetic moment around each 
Fe ion is found to be about 3.12/i#, which is well consis- 
tence with experimental results 3.33/i# at low tempera- 
ture 12K 19 . Furthermore, the calculated band structure 
of CAF state is shown that 2212* is a semiconductor with 
an energy band gap of around 0.64eV, as shown in Fig. 
[H(a), which is also well consistence with reported exper- 
imental results 19 . For 123*, the BAF state is the lowest 
energy state. In the BAF state of 123*, for each two-leg 
ladder, four spins group together become a super unit. 
Spins between two NN units are AFM. The coupling be- 
tween two NN two- leg ladders is also AFM. The moment 
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a). AFM 




d). BAF 



Chain Direction 

FIG. 3: (Color online) Schematic top view of five possible 
magnetic orders in the Fe-Fe square layer of the iron selenide 
compounds: (a) AFM Neel order in which the nearest neigh- 
boring Fe moments are AFM ordered; (b) Bicollinear-AFM 
(BCAF) order (the chain direction is changed into vertical 
direction for 2212*); (c) CAF order in which the Fe moments 
are FM ordered along the chain direction and AFM ordered 
across the chains direction; (d) BAF consisting of FM Fe4 pla- 
quettes tiled AFM along the chain direction; (e) Staggered- 
BAF (SBAF) configuration with FM diagonal double stripes 
that are also tiled AFM. 



TABLE II: Energetic and magnetic properties of the 2212* 
and 123*. Results in the magnetic states configurations, as 
shown in Fig. [3] using experimental crystal structure o 19 ^ 21 . 
AE is the total energy difference per iron atom in reference 
to the FM state, and rriFe is the local magnetic moment on 
Fe. 



2212* 



AE(eV)/m Fe (fi B ) 



FM 

AFM 

CAF 

BCAF 

BAF 

SBAF 



0/3.13 
0.2184/3.11 
-0.0118/3.12 

0.0944/3.13 



0.5 



-0.5 



123* 



AE(eV)/m Fe (fi B ) 



0/2.58 
-0.1131/2.38 
-0.1560/2.77 
-0.0139/2.55 
-0.1615/2.85 
-0.1514/2.74 
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around each Fe is about 2.85/i# for 123*, and the elec- 
tronic band structure calculated shows a semiconductor 
with an energy band gap £^=0.24eV, as shown in Fig. 
IH[b). The very small energy difference between FM state 
and CAF state is indicative of weak AFM coupling be- 
tween different chains, and the energy difference between 
AFM state and CAF state is indicative of strong FM 
coupling along the chain direction in 2212*. Similarly, 
in 123* system, the small energy difference between FM 
state and BCAF state is also indicative of weak AFM 
coupling between different ladders, and the energy dif- 
ference between BAF state and BCAF state indicative of 
four Fe atom plackets along the ladders. 

Nevertheless, these results are consistent with mag- 
netic exchange couplings obtained in other iron- 
chalcogenides FeTe and K0.sFe1.6Se2, where an FM NN 
coupling Ji, an AFM next nearest neighbor (NNN) cou- 
pling J2, and an third NN AFM J3 are necessary in de- 
scribing magnetic orders. As we will show in next sub- 
section, the magnetic orders of both materials can be 
obtained within models with the similar exchange cou- 
pling parameters. The AFM couplings between two NN 
chains in the 2212* and between two NN two- leg ladders 
in the 123* are exactly the third NN AF coupling, J 3 . 
The BAF order within each two-leg ladder can also be 
naturally understood from these couplings. Therefore, 
overall, the magnetism of all iron-chalcogenides can be 
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FIG. 4: (Color online) (a) Electronic band structure of the 
CAF state in 2212* with an energy band gap E g = 0.64eV; 
(b) Electronic band structure of the BAF state in 123* with 
an energy band gap .5^=0. 24eV. The Fermi energy is set to 
zero. 



unified into an effective model that includes local mag- 
netic exchange couplings as suggested in |32| and |33|. The 
values of exchange couplings can not be accurately de- 
termined since the result depends on the selection of the 
magnetic configurations^. 
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B. Magnetic Model for Ce 2 2 FeSe 2 

Following above results and effective models derived 
for other iron-chalcogenides, we construct the following 
model to describe this material 

# = JiJ^.^ + Ja Si-Sj + J'^Si'Sji (!) 

(ho) ((ho)) ihj) 

where J\ and J3 are the NN and next NN intrachain 
magnetic exchange couplings and J3 is the NN interchain 
magnetic exchange coupling, as shown in Fig. [3jc). 

The classical ground state of the Hamiltonian can be 
obtained exactly. In general, the classical energy is given 
as (for simplicity, we take S — 1) 

E c = '1J \ cos Q x + 4J 3 (cos 2 Q x - 1/2) + 2J' 3 COS U„, 

where (Q X1 Q y ) are the magnetic order wavevectors 
which can be viewed as the relative polarization angles 
between two NN intrachain spins and interchain spins re- 
spectively. The CAF phase is obtained when J3 < | Ji|/4 
and J3 > 0. 

In this state, we perform a linear spin wave analysis for 
this material in the classical limit. To do this we use the 
usual linearized Holstein-Primakoff transformation from 
spin operators to magnon operators which read as 

Sf = yf (bi + &!); SV = -iyf (&< " b\); S? = S - b% 

where i runs over all the lattice sites. Performing a fourier 
transform, the spin wave excitations of the model is given 
by 

k 

where, H = 2N J x + 2NJ 2 - 2NJ 3 is the ground state 
energy and S&\ = {b\,b_ k ), where 

&k — Ji(cos k x — 1) + J3(cos 2k x — 1) + J3 
bk = J3 cos k y 

Using the Bogliubov transformation, the linear spin- 
wave approximation, Eq. (|2]) can be diagonalized and 
shown in Fig. [5] It is interesting to see the effect of J3 
on the spin wave excitations. For J3 being AFM and close 
to 0.25Ji, the spin wave dispersion along chain direction 
(k X j7r) becomes quadratic at k x = 0. Otherwise, the 
dispersion is linear. 

C. Magnetic model for BaFe 2 Se3 

For the 123* system, we can begin with the following 
general Hamiltonian, 

H = J\b ^ Si ■ Sj + J\ a ^ Si ■ Sj + J[ a ^ Si • Sj 

+ J2 ^ Si ■ Sj + J3 ^ Si ■ Sj + J3 ^ Si • Sj, 

((id)) «<i,i»> <*,i> 




FIG. 5: (Color online) The spin wave dispersion relation 
as a function of k x -k y in the commensurate phase (0, ty) 
CAF state for 2212*. Parameters chosen are (Ji, J3, J3) = 
(-1,0.25,0.25). 

where Jib along with Ji a ,J( a and J3 denote the intral- 
adder vertical and horizontal NN and the third NN cou- 
plings, J 2 is the intraladder diagonal coupling and J3 is 
the interladder interaction as shown in Fig. [3jc). These 
coupling parameters reflect the symmetry breaking of the 
BAF state. 

We can treat the above model classically to obtain the 
exact ground state and phase diagram. We define the rel- 
ative polarization angles (Q x ,Q x ,Qy,Q y ) along the dif- 
ferent directions, with the primed variables going with 
the respective primed couplings. We can then write off 
the classical ground state energy as (for simplicity, we 
also take 5 = 1) 

E c = 2 J% cos Q y + 2 J lb cos Q y + 2J la cosQ x 
+ 2J[ a cosQ' x + 2J 2 cos Q y cosQ' x 
+ 2J 2 cos Q y cos Q x +4J 3 cos(Qx + Q' x ) 

We can then obtain the ground states by simply min- 
imizing the classical energy. With the BAF state being 
the ground state, we have (Q x , Q' x ,Q yi Q' y ) = (0, 7r, 0, 7r). 
Following the exchange coupling values measured for 
FeTe^ and K0.8Fe1.6Se2 35 , we expect that J\ a ~ Jib < 0, 
J 3 > 0, J3 > 0, J 2 > 0, and J' la > 0. The strength of 
the couplings satisfies, |Ji a | > J 2 > J3, J3, J{ a , which 
stabilizes the BAF phase. 

In the BAF state, we can obtain the spin wave excita- 
tions as done previously, which is given by 

where H = NJ la + NJ lb - NJ' la - 2NJ 3 - 
NJ' Z is the ground state energy and ^\ = 



5 




3.14 



-3.14 -3.14 



FIG. 6: (Color online) The spin wave dispersion relation 
of the lowest three branches (the other one is too high to 
be drawn in the same plot) as a function of k x -k y in the 
BAF state for 123*. The chosen parameters are fixed as 
(Jia, Jia, Jib, J 2 , J 3 , J'z) = (-1,0.1,-1,0.5,0.3,0.3). 



( b lk> 6 L b tk> b lk> K-ki 62,-fc, b Si -k, h-k)- 
Bk are four-by-four matrices, defined by: 



A k = 



I Eq Ji a J2 Jib 

Eq Jib J2 

Eo Jia 

\ . . . E 



Ak and 



(4a) 



2J 3 cos k x . The lower triangle elements are suppressed 
because both matrices are hermitian. 

By diagonalizing this Hamiltonian Eq. [3] for each k 
in the Bosonic metric, we obtain the spin wave disper- 
sion shown in Fig. [6] by taking Ji a = Jib = — 1, J2 = 
0.5, J3 = J3 = 0.3, J'i a = 0.1. The spin wave has four 
branches which is very similar to the BAF state discussed 
for KFe 2 Se2^. 



III. CONCLUSION 

In this paper, we have performed the first-principles 
calculations for the electronic band structures and mag- 
netic orders in these iron selenide systems featured with 
quasi-lD Fe chain (ladder) structures including 2212* 
and 123*. However, the calculated FS topologies still ex- 
hibit two dimensional or even three dimensional features. 
For 2212*, we find that the ground state is a CAF ordered 
semiconductor with an energy gap of 0.64eV, in agree- 
ment with the experimental measurements. For 123*, the 
calculated results show that the ground state is a BAF or- 
dered semiconductor with an energy gap of 0.24eV. These 
results suggest that that all magnetic structures in iron 
selenide compounds can be unified into an extended J\- 
J2-J3 model. We also calculate spin- wave excitations of 
the model which can be tested in future experiments on 
these two systems. 
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